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Abstract 

Q 

QQ In this article, we introduce iterative deterministic equivalents as a novel technique for the perfor- 

" mance analysis of communication systems whose channels are modeled by complex combinations of 

independent random matrices. This technique extends the deterministic equivalent approach for the study 

of functionals of large random matrices to a broader class of random matrix models which naturally 

arise as channel models in wireless communications. We present two specific applications: First, we 

^ consider a multi-hop amplify-and-forward (AF) MIMO relay channel with noise at each stage and 

^ derive deterministic approximations of the mutual information after the Kth hop. Second, we study a 

MIMO multiple access channel (MAC) where the channel between each transmitter and the receiver 

■«— j- is represented by the double-scattering channel model. We provide deterministic approximations of 

CN the mutual information, the signal-to-interference-plus-noise ratio (SINR) and sum-rate with minimum- 

mean-square-error (MMSE) detection and derive the asymptotically optimal precoding matrices. In both 

• scenarios, the approximations can be computed by simple and provably converging fixed-point algorithms 

> 

' ^ and are shown to be almost surely tight in the limit when the number of antennas at each node grows 

X 

infinitely large. Simulations suggest that the approximations are accurate for realistic system dimensions. 



The technique of iterative deterministic equivalents can be easily extended to other channel models of 
interest and is, therefore, also a new contribution to the field of random matrix theory. 
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I. Introduction 

Since the pioneering work of Tse and Hanly |^ on the capacity of code division multiple access 
(CDMA) technologies assuming long spreading sequences, the theory of large dimensional random 
matrices (RMT) has drawn an increasing interest from researchers in wireless communications and related 
fields lis, lO. RMT is in particular convenient for the study of multiple-input multiple-output (MIMO) 
channels H, 0, (random) linear precoders HI, ||6l, Q, multi-user systems 111,1191, multi-cellular systems 
ifTOl . lim . |[T2l . etc. In the early contributions, it was systematically assumed that the dimensions of the 
system under study could grow infinitely large and that the system performance admitted a deterministic 
limit that RMT can provide lUl, llH, |[T3l . It then became rapidly clear that, for most systems of practical 
interest, either the former assumption is not natural or the latter condition is not met. However, even for 
systems of large but finite size, the inherently random performance (e.g. instantaneous mutual information, 
signal-to-interference-plus-noise ratio (SINR)), can often be well approximated by deterministic quantities. 
Such quantities are called deterministic equivalents, and can be derived by various techniques, such as 
the Stieltjes transform method Q, Oil, the Gaussian method IHl, Ull, or the repUca method IHl, ifTSll . 

Deterministic equivalents are convenient to study the performance of wireless communication systems 
when a single system parameter can be modeled by a random matrix, e.g. the fading channel or the 
spreading codes. In order to tackle the performance analysis of more complex systems, such as multi- 
hop communications, random beamforming over random fading channels, or double-scattering channels, 
it is necessary to extend the notion of deterministic equivalents. In the present article, which is inspired 
by the original idea of |7|, where the performance of random isometric precoders over random fading 
channels is analyzed, we develop a systematic approach to generalize deterministic equivalents to iterative 
deterministic equivalents. To this end, we introduce a generic definition of deterministic equivalents of 
functionals of random matrices, which we extend, based on the Fubini theorem |fT9l . to a definition of 
iterative deterministic equivalents. 

As application examples, we then provide deterministic equivalents of the mutual information of the 
multi-hop amplify-and-forward (AF) MIMO relay channel ||20l, lUB, lUll (Section |III-A| ) and of the 
ergodic capacity as well as the sum-rate with minimum-mean-square-error (MMSE) detection of double- 
scattering multiple access channels (MACs) |[23l, | |24l (Section III-B| ). An overview of related research 
to both topics is provided in the respective sections. Our analysis is based on the Stieltjes transform 
method, documented in detail in |l3|- 

The remainder of this article is structured as follows. In Section |ll) we recall the fundamentals 
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of deterministic equivalents in RMT and develop the notion of iterative deterministic equivalents. In 
Section [Illj we study applications of iterative deterministic equivalents to the performance analysis of 
multi-hop relay channels and double-scattering MACs. The paper is concluded with Section |IV] All 
proofs, related results, and some exemplary Matlab codes are provided in the appendices. 

Notations: Boldface lower and upper case symbols represent vectors and matrices, respectively. 1m 
is the size-A^ identity matrix and diag(a;i, . . . , x^r) is a diagonal matrix with elements Xi. The trace, 
transpose, and Hermitian transpose operators are denoted by tr(-), (•)^, and (•)'^, respectively. The 
spectral norm of a matrix A is denoted by ||A||, and, for two matrices A and B, the notation A ^ B 
means that A — B is positive-definite. For a vector x = [xi . . . x^]^ , x > denotes Xj > for all i. 
The notations =^ and denote weak and almost sure convergence, respectively. We use CM (m, R) 
to denote the circular symmetric complex Gaussian distribution with mean m and covariance matrix R. 
We denote by the set [0,oo), by M" the set (— oo,0], and by i = -v/^. 1a{x) is the indicator 
function, i.e., 1a{x) = 1 iff x S ^ and 1a{x) = otherwise. E [•] denotes the expectation operator. For 
(«n)n>i and {hn)n>i two scqucnccs of random variables, we denote a„ x hn the equivalence relation 
o,n — bn -— > for n — >• oo. 

II. Iterative deterministic equivalents 

In this section, we will first recall the notion of deterministic equivalents in probability theory before we 
explain their connections to RMT and the performance analysis of communication systems. We will then 
introduce the Fubini theorem, which is the key ingredient to extend classical deterministic equivalents to 
iterative deterministic equivalents. 

A. Deterministic equivalents and random matrices 

Definition 1: Consider the probability space (17, F, P). Let (/n)n>i be a series of measurable complex- 
valued functions, /„ : 17 x C — )• C, and let {gn)n>i be a series of complex-valued functions, 5(„ : C — )■ C. 
Then {gn)n>i is a deterministic equivalent of (/n)n>i on Z) C C, if there exists a set yl C O with 
P{A) = 1, such that 

fniuJ,z) - gn{z) > 

n— >oo 

for all cj S A and for all z ^ D. 

Otherwise stated, a deterministic equivalent for (/n)n>i is a series {gn)n>i such that gn{z) approx- 
imates fn{uj,z) arbitrarily closely as n grows, for every z £ D and almost every uj. In particular, if 
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(/n)n>i converges almost surely to a limiting function /, i.e., for all {co,z) ^ Ax D with A C 
P{A) = 1 and Z) C C, we have 

fn{uj,z) >f{z) (1) 

n— >oo 

then {gn)n>i defined by = f, for all n, is also a deterministic equivalent of (/n)n>i- In many cases, 

one can further show that L /n(w, z)dP(Lo) — Qniz) > 0. Thus g„ is also an approximation of the 

expected value of /„. 

In the context of large dimensional random matrix theory, one often considers random matrices E 
i^Nxn growing dimensions A^, n — )• oo, where in general N/n = Cn is such that 

< liminf c„ < limsupc„ < oo. (2) 

" n 

This simply states that c„ is bounded so that the ratio N/n of the matrix dimensions is never too close 
to zero or infinity. Formally, to be in line with Definition [T] we will define random matrices in the 
following as series (H„)„>i = (H„(tj))„>i of matrices with growing dimensions which are defined on 
a probability space {Q,T,P), where every lo £ Q, generates the whole sequence (H„(a;))„>i and not 
only a single matrix H„(a;). 

In wireless communications, one is often interested in the behavior of functionals /„(H„,z), where 
H„ G i^Nxn ^ matrix describing the input-output relation of a wireless channel. In particular, 
fn{iln,z) = -^logdet(lAr + 2;H„HJ^), z £ R"^, is the (normalized) mutual information of the MIMO 
channel Jin between an n-antenna transmitter and an A^-antenna receiver at signal-to-noise ratio (SNR) 
z. Other quantities of interest are the SINR with linear detectors or precoders and the associated rates. 
The goal of a large system analysis based on RMT is to provide deterministic approximations of these 
random quantities, which become arbitrarily tight as the system dimensions grow. Thus, deterministic 
equivalents provide a deterministic abstraction of the physical layer. This is particularly interesting for 
complex channel models which are intractable by exact analysis. 

Deterministic equivalents for functionals of large dimensional random matrices have been considered 
for a wide range of communication channel models. For instance, in |5|, a deterministic equivalent for 
the ergodic mutual information of the Rician fading channel model = X„ + A„ is provided, where 
G d^A^xn j^^g independent entries with zero mean and a variance profile E[|(X„)jjp] = a'^^j, and 
A„ G (^Nxn ^ deterministic matrix. In ||25l . the deterministic equivalent of [5] is used to determine 
an asymptotically tight approximation of the ergodic capacity achieving input covariance matrix for the 
MIMO Rician fading channel. Deterministic equivalents were then extended to broader classes of wireless 
channel models, such as the capacity of the frequency-selective MIMO channel L15J . the MIMO MAC 
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with Kronecker correlation 11141 and the sum rate capacity of linearly precoded broadcast channels under 
imperfect channel state information Q- The application of such techniques is therefore very broad as it 
can simplify the difficult study of communication channels with a various number of random parameters 
(random channels, unitary precoders, path loss, etc.). Moreover, deterministic equivalents can be used to 
compute approximate solutions of otherwise intractable optimization problems ||T2|| . ||9l , ||25]| . 

All of the works mentioned above consider deterministic equivalents for random matrix models created 
from sums of independent random matrices. In many cases of practical interest, it is however necessary 
to consider more complex combinations of matrices, such as products or sums of products. These include 



for example the multi-hop relay channel (Section III-A i as well as the double-scattering channel model 



11231 (Section III-Bi. Another recent example is [7|, which considers random beamforming over fading 
channels, i.e., both the precoding and the channel matrices are assumed to be random. In this work, 
the authors derive deterministic equivalents of the mutual information and of the SINR with MMSE 
detection with respect to the random precoding matrices for quasi-static channels. Then, a second set of 
deterministic equivalents is found, treating both precoders and channel matrices as random. This technique 
relies on a fundamental result of probability theory, the Fubini theorem. In this article, we explain this 
approach in detail and generalize it to the new notion of iterative deterministic equivalents. 



B. The Fubini theorem 

Theorem 1 (SSS): Let (il, J", P) and P') be two probability spaces. Denote {Qx^', TxT , Q) 



their product space. Let / : x — )■ M be (J^ x J^')-integrable. Then 



f{uj,uj')dQ(uj,uj') = 



f{uj,J)dP'{oj' 
f{u,J)dP{u) 



dP{uj) 
dP'{J). 



In particular, consider a s&t A ^ F x F' . Then, we have from Theorem [T] that 



Q{A) 



uj')dQ{(jj, bj') 



lA{oj,bj')dP{uj) 



dP'U) 



(3) 



Equation Q is the core ingredient for the definition of iterative deterministic equivalents: Let (H„((^))„>i 
and (H^(a;'))„>i be two series of random matrices generated by the spaces {Q.,F,P) and {Q' ,F' ,P'), 
respectively. As in Theorem [T| call Q the product-space measure. Let /„((H„(a;), H'^(a;')), z) be a 



6 



functional of the matrices H„(a;) and H^(a;'). Assume that there is a function ^„(H„(a;), z), such that, 
for each uj e A C 0, with P{A) = 1, there exists a subset B{u}) C 17' with P'{B{uo)) = 1, on which 

/„((H„(a;),H;(a;')),z) -5n(H„H,z) ^0. (4) 

Although ^„(H„(w), is a random function (as it depends on w), it is independent of H'„(a;'). Thus, we 
can see ^„(H„(a;), z) as a deterministic equivalent of /„((H„(a;), H'„(a;')), z) with respect to (H'„(6t;'))„>i. 
Now, let us assume that there is a second function gn{z), such that f or uj £ C C 0, with P(C) = 1, 

gn(iln{uj),z) - gn{z) 0. (5) 

CdiWD = {{uj.uj') : uj £ AnC , J £ B{uj)} C 17 x il', the space on which /„((H„, H^), z)-5r„(z) 0. 
Then, from Q, this space has probability 

Q{D) =11 lD{^,u:')dP'{uj') 



(a) 
> 



(J 



Anc 



dP{oo) 

lD{uj,Lo')dP'{uj') dP{L0) 

B{u) 

dP{uj) 

Anc 

1 (6) 



where (a) is due to An B C ^ and B{oj) C 17', (6) follows since P'{B{uj)) = I iox uj £ A and (c) 
holds since n C) = P[A) + P(C) - P(A U C) = 1. 

To summarize, if a deterministic equivalent Qn exists for a functional of a random series (H^)„>i 
and a deterministic series (H„)„>i of matrices, and if additionally it can be proved that this deterministic 
equivalent holds true for almost every such (H„)„>i generated by a space 17, then the latter is also a 
deterministic equivalent for the random series ((H„, H'^))„>i. 

This is the mathematical key idea behind our method to derive iterative deterministic equivalents of 
functionals /„((Hn(a;), H'^(u;')), z), of two (or more) random matrices. First, one considers one of the 
sequences of random matrices, e.g. (H„(cj))„>i, to be deterministic and derives a deterministic equivalent 
with respect to (H'„(a;'))„>i. In the example above, this was the role of the functional 5„(H„(a;), z) 
which is independent of H^(cl!'). In a second step, one assumes the matrices (H„(cj))n>i to be random 
and derives an iterative deterministic equivalent gn{z) of g„(H„(Li;), z). Of course, this procedure can be 
carried out for any finite number of random matrices where in each step the "randomness" related to one 
of the matrices is removed. From the above construction, we will call {gn)n>i an iterative deterministic 
equivalent. 
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Fig. 1. Multi-hop amplify-and-forward MIMO relay channel. 

In the next section, we present two specific examples of iterative deterministic equivalents with 
applications to the capacity of multi-hop MIMO relay and double-scattering channels. From now on, all 
matrices and vectors should be understood as sequences of matrices and vectors with growing dimensions. 
For notational convenience, we drop the index n, e.g. we write H instead of (H„)„>i. 

III. Applications 

A. Multi-hop relay channel 

Consider a multi-hop AF MIMO relay channel where a source node communicates via K — 1 relays 
with a destination node. There is no direct link between the source and the destination and each relay can 
only receive data from the preceding hop. This is for example the case if the nodes follow a time-division 
multiple-access (TDMA) protocol where only one node is transmitting at any given time and the path 
loss between relay k and A; — 2 is large. Thus each data symbol reaches the destination after K channel 
uses. The source and destination are respectively equipped with n and uk antennas while the kih relay 
has nfe antennas. The relays operate an AF-protocoI where each node simply transmits a scaled version 
of its received signal to the next hop. We will consider a large system limit where n, ni , . . . , uk grow 
infinitely large at the same speed. Define the following quantities: 

n 

ci = — 
ni 

Cfe = ^, k = 2,...,K. (7) 

nk 

The notation "n oo" must be understood from now on as n ^ oo, such that < liminf^Cfe < 
lim sup„ Cjfc < oo for all k. We denote G C"*" the received base-band signal vector at the kih hop, 



8 



given by 

yi 




yfc = \/afcHfe J^^^yfe_i + rifc, k = 2,...,K (8) 

nk-i 

where G (^nkXUk-i ^ standard complex Gaussian matri}([^ (let no = n), X ~ CM{0, In) is the 
channel input vector, ~ CAA(0,I„^) is a noise vector, ak is a path loss factor, and the parameter 
is chosen to normalize the transmit power of the A;th node according to its power budget > 0, i.e., 

« _ ^0 _ 

PO — 1^1- TTT — PO 

-trE xx^ 

n L J 

The expectation in the last equation is with respect to the transmit and noise vectors onlyj^ The channel 
matrices and path loss factors ak are assumed to be known to the relays and the destination. Since the 
received signal at each relay is corrupted by noise, the system suffers from noise accumulation. This is 
in addition to the linear rate loss related to the TDMA protocol. Thus, the capacity decreases quickly 
with the number of hops K. Note that our system model is different from existing works which consider 
either no noise |[26l . or noise only at the destination lf22l . An exception is ll27l . in which the authors 
consider a similar system model, but do not provide closed-form expressions of the asymptotic mutual 
information. Several other works deal with the asymptotic capacity of the dual-hop relay channel ||28]| . 
|[29l . Recently, an exact expression of the mutual information of the dual-hop channel for finite channel 
dimensions was derived in |[30l . Here, we will provide an explicit deterministic equivalent of the mutual 
information at each relay for the general model ([8]). 

Let us introduce the following, recursively defined matrices [Pf,~i)- 



XX 



In 



Ro =E 

R, = E [yfeyHj = + ^h/^Uk^k-i {/3k-2) H^, k = l,...,K (10) 

and the functional J'k (^x,l3j.^i), x > 0, which are defined as 

Jfc(x,/3fc_i) = — logdetfl„, +x^^^^HfcRfe_i(/3fc_2)H^V k = l,...,K (11) 

'a standard complex Gaussian matrix X has i.i.d. elements Xij ~ CA/'(0, 1). 

^Under a long-term power constraint, the expectation could be taken also with respect to the matrices H^. Asymptotically, 
both constraints are equivalent (see Lemma |Tj(. 



9 



where /3;j = [/3o,--- ,/3fc]. With these definitions, we can express the normalized mutual information 
lk{0k-i) between and x as 

^kiPk-i) = ^ {Jk{h(3k-i)-Jk{hP'k-i)) (12) 
where P'f^ = [0,/3i, • • • ,(3^]. Next, we demonstrate by a simple example that (12 1 holds. 



Example 1 (2-hop Relay-channel): The normalized mutual information X2(/9i) between x and the 
channel output after the second hop y2 is given as 

= ^logdet L + k + ^H^hh)" ^^^^I^H^HiKHrh') 
Kn2 \ \ ni J mn j 

- ' logdet ( I., + ( I. + ^U,U^) U^) - -i- logdet ( I., + ^^H^RH 



Kn2 \ ^1 V ^ J J V ^1 

log det I Ins + 



^^H2Ri(/3o)hh) - ^ logdet ( I„, + ^^H2Ri(0)hh) 

ni J Kn2 \ ni J 



Kn2 \ ^ ni 

= 1(^2(1, /3i)-:^2(l,/3;)). (13) 

In the following, we will derive deterministic equivalents Tk0k-i) of ^ki(3k-i)- It will turn out that 
the recursive definition of the matrices Hkif^k-i) in ( [TOt allows us to calculate iterative deterministic 
equivalents of the mutual information after each hop. This is achieved by treating the matrix R/;,_i(/3;j_2) 
as deterministic and deriving a deterministic equivalent of J'k{x, (3i^_i) with respect to the matrix H^. 
This process can be iterated for Rfc_2(/3,fc_3), Rfc_3(/3fc_4), . . . and Hfc_i, Hjt_2, . . . until the determin- 
istic matrix Rq is reached. Before we address this problem, we will derive deterministic equivalents (3^ 
of the power normalization factors /J^: 

Lemma 1 (Asymptotic power normalization): Let /3o = /3o = Po- Then, 

h^^h = r—^ , k = l,...,K-l. 

n^oo 1 + akPk-1 

Proof: Recall the definition of Pk = j_trR ' ^bere R^ = R,fc For /c > 1, we have 

— trRfc = 1 H ""^ tr HfcRfc_iH^ 

Uk nkUk-i 

(a) 1 V 1 ~ [_| ~ 

= 1 + akPk-1 — / ^ hfc .-Rfc-ihfc,? 

nfc ^ Uk-i 

(^) , ^ Pk-i 1 , „ 

X 1 + aA;— ^ trRfc_i 

^trRfc_infc_i 

= 1 + akPk-i (14) 
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where (a) is obtained by denoting hkj £ C"''"^ the jth row vector of Hfc and (6) is due to Lemma [2] 
and Lemma |6] in Appendix [A| and the definition of /3fc_i. By the continuous mapping theorem llBTl . we 
finally have 

13k = -r^ ^ ^ . (15) 

-^-trRfc n^oo 1 + akPk-i 



In the next theorem, we will build upon Lemma [T] and provide deterministic equivalents of Ji^{x,(3k). 

Theorem 2: For k G {1, • • • ,K}, let I3j._i = [/3o ■ ■ ■ h-i] > be a sequence of random vectors, 
indexed by n, and = \Po - ■ ■ Pk-i] > be deterministic, such that /3j > ;Sj for i = 0, . . . , A; — 1. 

n— s>oo 

Then, 



Jk (a;,/3fc_i) - Jk {x,f3k-i) — ^ 

n— s>oo 

where Jk{x,Pi^_i) is recursively defined for > 2 as 



xakl3k-i 7, \ , 1 , xakf3k-i 



Jk{x,Pk-i) = CkJk-i ^ 1 — 7 — ^ ^,Pk-l +Cfclog 1 + 



Cfe + xafc^fc_i + efc_i (x,/3fc_i) j y Cfc + efe_i (x, /3fc_i) 

A , efc_i (a;,/3fc_i)\ ek-i{x,fik-i) 

+ log 1 + — — ^ r- 

\ Ck J Cfe + efe_i (x,/3fc_i) 

and Ck [x,0i^) for /c > is given by Theorem [3] The initial value J'l (x,^o) is given in closed form: 

^ f n\ 1 xaipo \ ( eo{xJo)\ eo(xjo) 
Ji [X, Po) = ci log 1 + + log 1 + — ^ , _ / 

Proof: The proof is provided in Appendix IB] I 



Theorem |2j allows us to compute the quantities J'k{x, (Sj^^i) recursively for any desired relay node 
k. The values of ek-i [x,(3f,_ij, needed at each stage, can also be calculated in a recursive manner as 
shown in the next theorem. 

Theorem 3: For k G {1, . . . , K — 1}, let = [/Sq - ■ ■ > be a sequence of random vectors, 
indexed by n, and = [/^o " " " f^k] > be deterministic, such that /3j > Pi for i = 0, . . . ,k. Let 
X > and denote by nik ix,(3k) = ^tr (ok+iPk^^'H.k+i'R-k (/3fe-i) H^+i + ^In.+i) ■ Then, 

i^k {x,Pk) - fhk {x,Pk} > 
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where nik [x,Pi.) is recursively defined for > 1 as 

XCk+l 



Cfe+i + efc (x,/3fc) 

and efc (x, is given as the unique positive solution to the following fixed point equation 

efc {x,^f,) = Cfe+i (cfc+i + efc (x,^^)) 

Cfc+i (cfe+i + gfc (x,^^,))^ _ / xak+ih a \ 

The initial values ?fio(x,^o) and eo(x,/3o) are given in closed form: 

mo{x, /3o) = —5 + (1 - ci)x 

Ci+eo(x,/9o) as- 



_ , xai/?o(l -ci) + ci V(xai/3o(l - ci) +ci) + 4xai/3oc? 
eo (x, /3o) = + . 

Proof: The proof is provided in Appendix |C] ■ 

Remark 3.1: The quantity mfc(x,/3^) can be seen as the Stieltjes transform |[3l of the empirical 
spectral distribution (e.s.d.) of the matrix ak+ifik:^ilk+i^k{Pk-i)^k-\-i evaluated at — ^. One can 
further show that Theorem js] implies the weak convergence of the e.s.d. ak+if3k^'H.k+i^k{Pk-i)'^k+i 
to a distribution function, whose Stieltjes transform is given by rhk, for almost every Hi, . . . 

Applying Theorem |2] and Lemma [T] to ([12]) yields the following corollary which provides a deterministic 
equivalent of the mutual information Ik{(3k-i)- 

Corollary 1 (Asymptotic mutual information of the K-hop AF MIMO Relay channel): 

-XfcOfc^i) k = l,...,K 

where 

{0k^i) = ^ (jk{h0k-i)-Jk{h0'k-i: 



K 

with = [h--- ^fc-i], ^Li = [0 ^1 • • • j3k-i] as given by Lemma[l} and Jk{x, Pk-i) and Jk{x, (3'k^i) 
as given by Theorem [2] 

Remark 3.2: The values of Tk{^k-i) '^an be very easily numerically computed. We provide the Matlab 
code which was used to generate the numerical results in this section in Appendix |j] Due to the recursive 
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implementation, the computational complexity grows quickly with k. Calculating Jk{x, with high 

precision for large values of k (> 10) seems therefore impractical. 

We would now like to verify our analysis by some numerical results. To this end, we consider a system 
with three relays, i.e., K = A. We assume n = 71,4 = 4, ni = 713 = 8, n2 = 12, pi = ps = 0.7po and 
P2 = 0.5po- The last assumption allows us to control the transmit power of all nodes by the transmit 
SNR pq of the source node. We further assume the path loss factors ai = 1, 02 = = 0.7, 03 = 0.5. 
Fig. [2] shows the average normalized mutual information E [^2^fc [Pk-i]] ^f^i" each hop (/c = 1, . . . , 4) 
versus the transmit power po of the source node. Note that we have re-normalized all results by ^ to 
put them on a common ground for comparison. The deterministic equivalents {f^k-i) provided 
by Corollary [T] are drawn by solid lines, simulation results are represented by markers. The error bars 
represent one standard deviation of the simulation results in each direction. We can observe a very good 
fit between the asymptotic approximations and the simulation results for all k and the entire range of pQ. 
As expected, the performance decreases rapidly with each hop. 

Finally, we would like to remark that, although we have considered a rather simple channel model 
with neither antenna correlation nor precoding at the nodes, more involved channel models can be treated 
in a straightforward fashion with the same techniques. 

B. Double-scattering MAC 

Consider a discrete-time MIMO MAC from K transmitters, equipped with n^, k = I, . . . , K, antennas, 
respectively, to a receiver with antennas. The channel output vector y G reads 

K 

y = HfcXfc + n (16) 

k=l 

where S C^^"'' and xj. G C"*- are the channel matrix and the transmit vector associated with the 
kih transmitter, n ~ CAf{0, ^Ia?) is a noise vector and p > denotes the SNR. We assume Gaussian 
signaling, i.e., x/j = [s^fc^i, • . • )2;fc,nfc 

y ~ CJ\f{0,Qk), where Qfc G C"''''"\ The channel matrices 
are modeled by the double-scattering model ||23]| 

where G C^""^, Sk G C^*''^'' and G C"""""" are deterministic correlation matrices, while 
Wi^fc G C^^^*" and W2.A; G C^*^"*" are independent standard complex Gaussian matrices. Since the 
distributions of Wi ^ and 'W2^k are unitarily invariant, we can assume = diag(sfc,i, . . . , s^^Nk) to be 
diagonal matrices, without loss of generality for the statistics of y. 
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Fig. 2. Average normalized mutual information E [^Xfc(/3j._j^)] after the fcth hop versus the transmit SNR po of the source 
node. The deterministic equivalents ^Xfc drawn by solid lines, the simulation results by markers. The error bars 

correspond to one standard deviation of the simulation results in each direction. 



The double-scattering model ||23l was motivated by the observation of low-rank channel matrices, 
despite low antenna correlation at the transmitter and receiver, see e.g. ll24l . ll32ll . A special case of the 
double-scattering model is the keyhole channel ll33l . ll34l . which exhibits null antenna correlation, i.e., 
Rfc = In and = for all k, but only a single degree of freedom. The existence of such channels 
(under laboratory conditions) was confirmed by measurements in ll34l . Several theoretical works have 
studied the double-scattering model so far. The authors of |35] derive capacity upper-bounds for the 
general model and a closed-form expression for the keyhole channel. An asymptotic study of the outage 
capacity of the multi-keyhole channel was presented in |l36l. The diversity order of the double-scattering 
model was considered in ll37l and it was shown that a MIMO system with t transmit antennas, r receive 
antennas and s scatterers achieves the diversity of order trs/ max(t, r,s). A closed- from expression of the 
diversity-multiplexing trade-off (DMT) was derived in ll38l . Beamforming along the strongest eigenmode 
over Rayleigh product MIMO channels, i.e., the double-scattering model without any form of correlation, 
was considered in f39l. Here, the authors derive exact expressions of the cumulative distribution function 
(cdf) and the probability density function (pdf) of the largest eigenvalue of the Gramian of the channel 
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matrix and compute closed-form results for the ergodic capacity, outage probability and SINR distribution. 
In a later paper BOl . the MIMO MAC with double-scattering fading is analyzed. The authors obtain 
closed-form upper-bounds on the sum-capacity and prove that the transmitters should send their signals 
along the eigenvectors of the transmit correlation matrices in order to achieve capacity. Despite the 
significant interest in the double-scattering channel model, little work has been done to study its asymptotic 
performance when the channel dimensions grow large. We are only aware of ll32ll . in which a model 
without transmit and receive correlation is studied relying on tools from free probability theory. Implicit 
expressions of the asymptotic mutual information and the SINR with MMSE detection are found therein. 

In the following, we provide deterministic equivalents of the mutual information, the SINR with MMSE- 
detection and the associated sum-rate. In addition, we derive the precoders which maximize the determin- 
istic equivalent of the mutual information and provide a simple algorithm for their computation. The key 
idea behind the following proofs is that the double-scattering channel can be interpreted as a Kronecker 
channel |[T4l with a random receive correlation matrix, which itself is modeled by the Kronecker model. 
This observation allows us to build upon |[T4l which provides an asymptotic analysis of the performance 
of Kronecker channels with deterministic correlation matrices (Theorem |9] in Appendix |A]|. Based on 
the Fubini theorem, we extend this work by allowing the correlation matrices to be random. A similar 
technique can be applied to more involved channel models, such as channels with line-of-sight components 
or MIMO product channels with an arbitrary number of matrices. 

Denote In{p) the instantaneous normalized mutual information of the channel (16 1, defined as iHTI 

In [p] = ^ logdet ^I^ + p J^HfcQfcH^^ . (18) 

Moreover, denote ■ (p) the SINR at the output of the MMSE detector related to the transmit symbol 
Xk,j, given by EH 

K 



{p) = ( ^ H,hH - + -In I h^,,. (19) 



We define the normalized sum-rate Rn {p) with MMSE detection as 

= ]^ E E (1 + <j ^p)) ■ (20) 

k=l j=l 

The notation 'W — )• oo" will be used to denote that N and all Nk, grow infinitely large, satisfying 
< liminf ^ < limsup=^ < oo and < liminf ^ < limsup^ < oo. Additionally, we need the 
following technical assumption: 
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A 1: For all k, limsup^||Rfc|| < oo, limsup^||Sfc|| < oo and limsup^||TyfcQfc|| < oo. 

Remark 3.3: This assumption implies in particular that the antenna correlation at the transmitter and 
receiver side cannot grow with the system size, as it would be the case for very dense antenna arrays 
ll43l . Amendments to relax this assumption can be made, following the work in |[T4l . Moreover, the last 
constraint, limsup;v||TjtQfc|| < oo, implies that no transmitter is allowed to focus an increasing amount 
of transmit power in a single direction. 

Our first theorem introduces a set of 3K implicit equations which uniquely determines some quantities 
{gk,gk,^k) (1 < A; < K). These will be needed later on to provide deterministic equivalents of In{p), 
7fc^,.(a2),and Rn (p). 

Theorem 4 (Fundamental equations): The following system of 'iK implicit equations in g^, Qk, and 

5k{l<k< K): 

5fc = -^trT|Q,T|(5feT|QfcT|+I„,^ ' 



1 X 

9k = —2_^ r—^ ^ (21) 

s -1 



has a unique solution satisfying gk,gk,^k > for all k and p > 0. 

Proof: The proof is provided in Appendix |D] ■ 

Remark 3.4: The values of g^, gk, and 5k can be computed by a standard fixed-point algorithm which 
iteratively computes ( [2T] ), starting from some arbitrary initialization 9k'\gk'\^k'^ > 0- This algorithm is 



proved to converge, generally terminates within a few iterations (depending on the system size and the 
desired accuracy), and does not pose any computational challenge. 

The next theorem provides a deterministic equivalent of the (ergodic) mutual information based on the 
quantities {gk,gk,Sk) as provided by Theorem |4] 

Theorem 5 (Mutual information): Assume that A 1 holds. Then, 

(i) /^(p)-I^(p)-^0 

N-^oo 

(ii) E[In(.p)]-In(.p)- ^0 
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where 



K 



k=l 



1 ^ 

- ^ l^logdet (Iat, + gkSkSk) + logdet (in^ + 5fcT^QfcT|^ - 2nkgkgk 



k=l 



(22) 



and Qk^Qki ^k are the unique positive solutions to ([21 

Proof: The proof is provided in Appendix [E] ■ 

The following result allows us to compute the asymptotically optimal precoding matrices which 
maximize In {p) under individual transmit power constraints. 

Theorem 6 (Optimal power allocation): The solution to the following optimization problem: 

(Q^ , . . . , Q^) = arg max In {p) 

s.t. — tr Qfc < Pk yk 

nk 

where Ijy (p) is defined in Theorem 5 is given as = UfcP^U^, where Ufc G C^'^^^^" is defined by 
the spectral decomposition of = Ufcdiag(tfc^i, . . . ,ffc^„jU^ and = diag(p^ ^ . . . ,p^ „J is given 
by the water-filling solution: 

Ph = f/^fc - T^^l (23) 



\ 9ktk,j / 

where jjk is chosen to satisfy ^trP^ = and g'^ is given by Theorem |4] for = Q^. 

Proof: The proof is provided in Appendix |F] ■ 

Remark 3.5: The optimal power allocation matrices P^ can be calculated by the iterative water-filling 
Algorithm 1 (see llT4l Remark 2] and Q Remark 3] for a discussion on the convergence of this algorithm). 



Remark 3.6: Denote by (Q^, . . . , Q^) the set of precoding matrices which maximize E [In{p)] for a 
given set of power constraints. If the condition limsup||TyfcQ^|| < oo holds for all k, then E [In{p, Q*, • • • , Q 
77v(p, Qi, . . • , Q*k) ^ 0' by Theoremp^and the strict concavity of In{p) and In{p) in the matrices 

Af— 5>oo [_J 

Qfc. However, this condition is difficult to verify and is outside the scope of this paper. See [25] for such 
a technical discussion in the case of Rician fading channels. 

Next, we provide deterministic equivalents of the SINR (p) at the output of the MMSE detector 
and the associated sum-rate Rn (p)- 
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Algorithm 1 Iterative water-filling algorithm 



1: Let e > 0, n = and p^'^ = Pk for all A:,j. 
2: repeat 

3: For all k, compute (7^'" according to Theorem 



4: For all k,j, calculate p^'"^^ 



with matrices = U^diag (p^'") U^. 



.rt) ' such that ^ E -iiP^f' = 



n 



n+1 



6: until maxfcj |p 



•->t,n 



-*,n— li ^ 



Theorem 7 (SINK of the MMSE detector): Let = diag {pk,i, ■ ■ ■ ,Pk,nk) ^^id = diag(tfc^i, . . . , tfc,„^) 
for all k. Assume that A 1 holds. Then, 

where j^j (p) = Pkjtkjgk and is by given by Theorem |4[ 

Proof: The proof is provided in Appendix |G] ■ 

Remark 3. 7: It is easy to see that the theorem is also valid under the more general assumption = 
Ufcdiag(tfc,i, • • . ,tfc,njU^ and Qfc = Ufcdiag(pfc,i, . . . ,pfc,„jU^. 

Corollary 2 (Sum-rate with MMSE decoding): Let = diag {pk,i, ■ ■ ■ ,Pk,nk) and = diag(tfc^i, . . . , t^^nk] 
for all k. Assume that A 1 holds. Then, 

{ii) E[Rn{p)]-Rn{p) 

Af-s>oo 

where 

K nu 



N 

k=i j=i 



and the 7^^- (p) are given by Theorem jvj 

Proof: The proof is provided in Appendix IH] 



Remark 3.8: Careful inspection of i\22\ reveals that the third term of In (p) equals Rn{p) since 



- J]logdet(l„, + <7fcT|QfcT|) = - J] J] log (1 + p,,,*^,,^^) . (24) 



k=l k=l j=l 
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Thus, all other terms in p2] ) correspond to the gains of successive interference cancellation [,44.| over 
simple MMSE detection. 



A special case of the double-scattering channel is the Rayleigh product MIMO channel 113911 which 
does not exhibit any form of correlation between the transmit and receive antennas or the scatterers. For 
this model, Theorems [4j [5] and [7] can be given in closed form as shown in the next corollary. 

Corollary 3 (Rayleigh product channel): For all k, let N/. = S, rij. = N and assume = I^v, 
Sfc = I5, R/c = liv, and = 1^. Then In (p) and (p) as defined in Theorems js] and jv] can be 
given in closed form as 

lNip) = log (^ + P^9 + f (^-1)) -Klogig)-2K{l-g) 



and 



where g is the unique root to 



7k,j (p) = — — 



s 2f S 1\ ( S 1 S ( l\\ SI 

such that 5 G (1 — min [;^, , l). 

Proof: The proof is provided in Appendix |l] ■ 

Note that similar expressions for the asymptotic mutual information and MMSE-SINR have been 
obtained in |[32l by means of free probability theory. However, these results require the numerical solution 
of a third order differential equation. 

As a first numerical example, we consider the "multi-keyhole channel", i.e., K = 1, Si = 
Ri = Iat, Ti = Qi = for = ni = 4. Fig. |3] depicts the normalized ergodic mutual information 
E [In{p)] and its asymptotic approximation In{p) versus SNR for different numbers of "keyholes" Ni G 
{1,2,3,4,100}. Surprisingly, the match between both results is almost perfect although the channel 
dimensions are very small. As one expects, the multiplexing gain increases linearly with Ni until Ni > 
N = 4. Larger values of A^i only change the statistical distribution of the channel matrix while the degrees 
of freedom are limited by the number of antennas (for A^i — )■ 00, Hi becomes a standard Rayleigh fading 
channel l|23l ). 

As a second example, we consider a MAC from K = 3 transmitters, assuming the double-scattering 
model in [23 1. Under this model, the correlation matrices are given as R^ = G(0r,/c, c^r.fc, A^A:)> Sfc = 
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-10 -5 5 10 15 20 25 

P (dB) 

Fig. 3. Ergodic mutual information E [/jv(p)] of tlie multi-keyhole channel and its deterministic equivalent In{p) versus p. 

G{4>s,k,ds,k,Nk) and = G{(t)t,k,dt,k, Nk), where G{(t),d,n) is defined as 

[Gicl>, d, ^ = ^ I] exp (^i27rdik - I) sin (y^) ) ' ^^^^ 

The values (pt^k and (pr^k determine tlie angular spread of the radiated and received signals, dt^k and 
dr^k are the antenna spacings at the kth transmitter and receiver in multiples of the signal wavelength, 
Nk can be seen as the number of scatterers and ds^k as the spacing of the scatterers. For simplicity, 
we assume N = A, = l/n^, = 11, Uk = 3, dt^k = c?r,fc = 0.25 and dg^k = 50 for all k. We 
further assume (f>j.^f^ = cpt^k for all k, with (pr^k G {^/4, vr/2, vr} and (f>s,k = "^/S- Fig. |4] shows E [I^ip)] 
and In{p) with uniform and optimal power allocation versus SNR. Again, our asymptotic results yield 
very tight approximations, even for small system dimensions. Note that we have used the precoding 
matrices provided by Theorem [6] for the simulations as the optimal precoding matrices are unknown. 
For comparison, we also provide the sum-rate with MMSE detection '&[Rn{p)\ and its deterministic 
approximation Rj\^{p). We observe a good fit between both results at low SNR values, but a slight 
mismatch for higher values. This is due to a slower convergence of the SINR ^y^jip) to its deterministic 
approximation J^jip)' ^^^^ documented in the RMT literature, e.g. P31 . 
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P (dB) 

Fig. 4. Ergodic mutual information E[/jv(p)] and sum-rate E[i?jv(p)] of the multiple access channel and their asymptotic 
approximations /jv(p) and Rn{p) versus p. 

IV. Conclusion 

In this paper, we have presented a novel tool for the large system analysis of communication systems, 
called iterative deterministic equivalents. This tool is particularly suited for the analysis of channel models 
composed of complex combinations of independent random matrices, e.g. products or sums of products 
of matrices. We have demonstrated the usefulness of this approach with the help of two examples which 
had not been solved in the literature before. These are the multi-hop AF MIMO relay channel with 
noise at each stage and the MIMO MAC under the double-scattering channel model. For these channel 
models, we have provided asymptotically tight deterministic approximations of information theoretic 
quantities of interest, such as the mutual information and the sum-rate with MMSE detection. These 
approximations can be easily computed by provably converging fixed-point algorithms and do not require 
any numerical integration. Our simulation results suggest that the asymptotic performance approximations 
are very accurate for finite system dimensions with only a few antennas at each node. Finally, the method 
of iterative deterministic equivalents is applicable to a wide range of channel models of interest (e.g. 
combinations of correlated i.i.d. and random unitary matrices Q) which cannot be easily treated so far 
with other techniques. 
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Appendix A 
Related results 

Lemma 2 ([46. Lemma 2.7], [9 Lemma 4]): Let (A^fj^^^ be a sequence of random N x N matrices, 
satisfying limsup^|| A^vH < oo, almost surely. Let (x7v)/yr>i be a sequence of random A^-dimensional 
vectors of i.i.d. entries with zero mean, variance and 8th order moment of order 0(l/A^^), 

independent of Atv- Then, 

H A ^ i * a.s. „ 

x^Aatxat - ^tr Aat > 0. 

Lemma 3 (Matrix inversion lemma M2 Eq. (2.2)]): Let A G c^""^ be Hermitian invertible. Then, 
for any vector x G and any scalar r G C such that A + rxx"^ is invertible, 

-1 x^A-i 



x^ (^A + rxx^ 



1 + TX^A^^X 



Lemma 4 (Rank-l perturbation lemma l{47\l): Let z < 0, A G C^^^, B G C^^^ with B Hermitian 
nonnegative definite, and v G C^. Then, 

lAII 



tr ( (B - zIn)~^ - (B + vv" - zIn)-^ ] A 



< 



Lemma 5: Let R G C^^^ be Hermitian with smallest eigenvalue Amin > 1 and a, b,c,d > 0. Then 
\trR(aR + 6lAr)"^R(cR + (fIiv)"^ > ^ 



Proof: Let R = UAU^, where U G C^""^ is unitary and A = diag {6i,...,6n)> 1- Thus, 
-^tr R (aR + 5Iiv)"^ R (cR + dlN)'^ = ^tr (aA + bl^v)"^ (cA + dl^v)"^ 

1 1 ^ 5Ua + b){c + d) 



{a + b){c + d)N {a6i + 6)(c(5i + d) 
- {a + b)ic + d)- ^^^^ 



Lemma 6: Let the matrices Rfc(/3;,_^), be defined as in (lOi. Then, almost surely 



limsup||Rfc(/9fc_i) II < oo, /c = 0, . . . , -R'. 
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Proof: For G {1, . . . , K}, denote by {Vik-, Fk^ Pk) the probability space generating the sequences 
of random matrices Hfc. By P6l . we have on a space Bk C ilfc with Pk{Bk) = 1, 



1 + 



0. 



Obviously, we have ||Ro| 



1. Thus, almost surely. 



limsup ||Ri(/3o)|| < 1 + limsup 



ai/3o 



n 



HiH^ 



1 + ai/3o limsup ( 1 + 



< CX). 



(28) 



(29) 



Consider now the product space (f^i x 0.2,^1 x J^2, Q2)- By the Fubini theorem, we have on a subspace 

C2 C f^l X ^2 with Q2{C2) = 1, 

limsup||R2(/3i)|| < 1 + limsup H2Ri(/3o)H^ 

1 



< 1 + limsupa2/3i||Ri(/3o)||- 



ni 



H2H2 



1 + a2/3i limsup 1 + ai/3o 1 H = 



1 + 



< 00. 



(30) 



Repeating the last step k — 2 times concludes the proof. 



Definition 2 (Standard interference function MH/j-' A function h(xi, . . . , xk) = [hi{xi^ . . . , x^), . , 
hxixi, ■ . . ,xk)]^ G is said to be standard if it fulfills the following conditions: 

1) Positivity: for each j, if xi, . . . , xk > 0, then hj{xi, . . . , xk) > 0. 

2) Monotonicity: if xi > x']^, . . . , Xi^ > x'^, then for all j, hj{xi, . . . , xk) > hj{x'i, . . . , x'^). 

3) Scalability: for all q > 1 and for all j, ahj{xi, . . . , Xi<:) > hj{axi, . . . , ax^c). 



Theorem 8 (Fixed-point theorem l[48\ Theorem 2]): If a ii'-variate function h(xi, . . . ,xk) is standard 
and there exists (xi, . . . , xk) such that for all j, xj > hj{xi, . . . ,xk), then the fixed-point algorithm 
that consists in setting 



for t > 1 and for any initial values x^''^ . . . , x^^ > converges to the unique jointly positive solution 
of the system of K equations 



Xj = hj{xi, . . .,xk), j G {1, • • ■,K}. 
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Theorem 9 (1(14] Corollary 1, Theorem 2]): For k G {\, . . . ,K}, let {nk)^^^ = {'nk{N))^y-y be 
a sequence of positive integers and let (Rfc^jv)^^^ R^^at G C^^^, {Tk,N) N>v'^k,N & C""=^"'\ 
and (D7v)7v>i> G C^^^, be three sequences of nonnegative definite Hermitian matrices, sat- 
isfying limsupjv||R.fc,Af II < c»> limsupjv||Tjt^jv|| < oo, and limsup;v||DAr|| < cxd. Let {^k,N) ^yi^ 
^k,N G C^^"*", be a sequence of random matrices with i.i.d. complex Gaussian entries with zero mean 
and variance l/n^. Denote Bat = J^k'^k N-^>',NT^k,N^k N'^k N define the function Vn{x) = 

logdet (liv + xBjv) for x > 0. Let Ck = rik/N and assume that < liminfjv Ck < limsup^r c^t < oo 
for all k. Then, 

(i) ^tr T>N (bn + -In) - ^tr Djv f V CiKi^N + -In] 



(u) Vn{x)-Vn{x)-^Q 



where 



1 f ^ \ ^ 1 1 ^ 

— logdet I Iat + X ^ ek,N'Rk,N 1 + XI ]v ^^^^^^ + ^k,NT^k,N) " ]^ XI ' 

V k=l / k=l k=l 



Vn (x) = — logdet ( Iat + X 2^ ek,N'Rk,N ) + X ^ logdet + ek,NT^k,N) - nkek,Nek,N 
and where ^k,N^ k = 1, . . . ,K, we given as the unique solution to the equations 

&k,N = trTfc TV (efc,ArTfc TV + In J ^ 

nk 

1 1 

ek,N = — trRfc.Tv y^eiRi^N + -Ia^ 
such that ek,N, ^k,N > for all k. 



Corollary 4 (Special case of Theorem^ see also [47]): Let (^^)Tv>l = i^{^))N>i ^ sequence of 
positive integers. Let (Rtv)tv>i' I^n £ C"^", be a sequence of nonnegative definite Hermitian matrices, 
satisfying limsupTv||RAf || < oo and let (Xtv)tv>1' ^ C^^", be a sequence of random matrices 
with i.i.d. complex Gaussian entries with zero mean and variance 1/n. For x > 0, define the following 
functions rriTv (x) = ^tr (XtvRatX^ + ^Iat) ^ and Jtv (x) = ^ logdet (Itv + xXatRtvX^). Denote 
c = ^ and assume that < liminfTv c < limsupTv c < oo. Then, 



where 



(i) mN{x) — mN{x) — ^— 0, {ii) Jn{x) — Jn{x) — > 



_ . . 1 / Rtv Ix 

mN{x) = — tr , _ + -In + (1 - c)x 
iV \c + eTv X J 

Jn{x) = ^ logdet ([c + cat] I„ + xRat) + (1 - c) log (c + cat) ^ log(c) 

iV c + eTv 
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and eiv is defined as the unique positive solution to the implicit equation 



1 „ / Rtv 1. ^ 



eN = ^trRN[—^ + -ln] ■ (31) 

A* \c + cn X 



Appendix B 
Proof of Theorem[21 



First, notice that 



< Oik 



nk-i 



HfcRfc„i(/3fc_2)H^ ^ HfcRfc_i(/3fc_2)Hfc 



nk-i 
HfcHfc 

nk~i 



rik-i 

|/3fc_iRfe_i(/3,,_2) - Pk-iTik-iiPk-2) 



k-l 



/3fe-l| Rfc-l(/3fe-2) + l^k-1 Rfc-l(/3A:-2) ~ Rfe-l(/3fc-2) 



(32) 



Since, almost surely, lim sup II Ri(/3o) — Ri(/3o)|| < limsupai|/3o- 



and limsup ||Rfc_i(/3^_2)|| < 



oo (see proof of Lemma [6]), one can iteratively show that r]k — >■ 0, almost surely. Thus, 



\Jk {x,f3k-i) - Jk {x,(3k- 



0. 



(33) 



This means that we can from now on replace by (3k and focus on Jk {x,(3k~i)- 

As a consequence of Corollary |4| Lemma |6j and the Fubini theorem, we obtain the following relation 



Jk{.x,~P, 



k-lj 



Jk{x,P, 



k-l) 



k > 1 



(34) 



where 



1 



Jk {x,Pk-i) = — logdet ([cfc + ek-i {x,Pk-i)] In^-i + xafc/3fc-iRfc-i {Pk~2)) 

IT'k 

+ (1 - Cfc) log {ck + Bk-i {x, fik-i)) - ^^/^^^^'^ ^ - log (cfc) (35) 



Cfc + efc_i [x,(3k-i) 



and efc_i {x,(3k-i) is given as the unique positive solution to 



ek-i {x,Pk-i) = — trafc/3fc_iRfe_i {(3k-2) 

^k 



Cfc + Ck-l (x,Pk-i) ^ 



(36) 



In particular, for /c = 1, we have 



Ji {xjo) =Jl {xjo) 



(37) 
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where 



Jx (x, ;3o) = ci log 1 + 



xai/3o 



ci + eo (x, /3o) 



+ log 



^ ep (x, /3o) \ _ ep (x, /3o) 



ci y ci + ep (x, /3p) 



_ , xai;So(l - ci) +ci a/ (xai/3p(l - ci) + ci) + 4xai/3pcf 
ep (x, /3p) = ^ V 



(38) 



(39) 



2 2 
according to Corollary |4] Note that ( |3T] ) permits a closed-form solution ( |39l ) in this case. 

Now, using the recursion Rfc_i = + ^^^^^i^Hfc_iRfc_2 (^^.g) H^,^ ([TOll in (|35]), we 

obtain 



Jk{x,l3k-i) = CkJk-i 



efc-i (x,^fc_i 



>/3fc-2 +Cfclog 1 + 



\ Ck J Cfc + efc_i (x,/3fc_i) 

In the proof of Theorem [3} it is shown that 

ek-i (x,^fc_i) - efc_i (x,^fc_i) 0. 

By the continuous mapping theorem IIBTTI . we therefore have 

^ ( xakh-i \ 

I 1 + xafe/3fc_i + efe_i (x,/3fe_i) ' ''^'^ J 



(40) 



(41) 



Jk-l 



xuklS, 



k-1 



-> 0. 



(42) 



Applying the last result together with Corollary |4} Lemma [6j the continuous mapping theorem and the 



Fubini theorem to ( |40| ) concludes the proof for /c = 2 since (x,/3o) = J^i (x,/3o) by ( |37] ). The 
convergence for A; > 2 is shown by successive iterations of the last steps. 



Appendix C 
Proof of Theorem[3] 



From standard matrix inequalities and (32i, it follows that 



\mk (x, Pf,) - ruk (x, I < x^ak 



+1 



nk 



I ARfc - h^k iPk-i) II 0. (43) 



Thus, we can replace from now on /3fc by /3k, for almost every (Hi, . . . , Hr)- 
From Corollary |4j Lemma [6] and the Fubini theorem, it follows that 



ITT-k {x,Pk) - l^k {x,(3k} 







(44) 
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where 

mk{x,Pk) = tr / / + -In + (1 - Ck+i)x (45) 

and efc (x, /Sj^.) is given as the unique positive solution to 

e, {x,-Pk) = — tra,^.i/3,R, (^>!±l^^^^i^ + J . (46) 

nk+i \Ck+i+ek{x,Pk) X J 

In particular, we have mo{x,0k) — ""^0(2;, /3fc)' where 

rhoixJo) = 3 ^ + (i_ci)x (47) 

Ci+eo(a;,,9o) ^ 



xai/3o(l -ci) + ci V(xai/3o(l - ci) +ci)^ + 4xai/3ocf 
eo (a;, /3o) = ^ h -'^ . (48) 

Replacing {^k-i) in ^ by its recursive definition {0k-i) = ^n. + ^^^^^Hk^k-i {Pk-2) ^k 
To] ) yields after straightforward calculus 



e/c {x,Pk) = Cfc+i (cfc+i + efc {x,(3k)) 

Cfc+i (cfc+i + efc / xofc+i^fc A 

^ =^ "mk-l = 7 ^^)Pfe-l • (49) 

xak+iPk yck+i+xak+if3k + ek[x,(3k) J 

Similarly, one obtains 

~ ^ z, ^ Ck+i{ck+i + ek{x,Pk)) I xak+ih a \ . l^ \ 

mk{x,(3k) = = ruk-i =— 7 — ^,Pk-i + (1 - Cfc+i)x. 

Q-k+iPk \ Cfc+i + xak+iPk + ek [x, /3k ) I 



(50) 



Combining the last two equations leads to 



fhkix,(3k) = ^^''^j - . . (51) 

Consider now the quantity (x,3a;)> k>\, defined as a positive solution to 

efc {x, Pk) = Cfc+i (cfc+i + ek (x, Pk)) 

Cfc+i (cfe+i +efc (x,^fc))^ _ / xak+ih a \ 

xak+\Pk y Cfc+i + xafc+i/3fc + ek (x, j 

where fhk (x, /fl^) is recursively defined for /c > 1 as 

^k {x,Pk) = ^^-^j -r, N • (53) 



27 



It remains to show that a unique solution to ([52]) exists and that ek (x,/3u) — Ck (x,Pj^) > 0. Let us 



first define the following functions for A; > 1: 

fk{z) = Ck+1 (cfe+1 + z) = nik-i ■ ^— — ,Pk-i (54) 

I r \ I , ^ Cfe+i (Cfc+i + zf _ ( Xttk+ih 7, \ 

Jk{z) = Cfc+i (cfc+i + z) = nik-i ■ ^— — ■ (55) 

xak+iPk \Ck+i + xak+iPk + z J 

From ( [46l ) and with the help of Lemma |5] (note that the smallest eigenvalue of is greater or equal 
to 1 for all k), one can easily verify that fkiz) satisfies the following properties for z > 0: 

(i) 

1 =- > 



Ck+i + z + xak+il3k 



fk{z) > Ck+l{Ck+l + z) 

(ii) for z> z' >{), 

fk{z) - fk{z) 

{z - z')ck+ial^Jl 1 f ak+iPkTik , 1^ f ak+iPkTik , 1 

— 7 — 7\7 ; — \ — trUfc ; 1 Irifc tlk ; — r ^ 

[Ck+l + z'){Ck+l + Z) Uk \Ck+l + Z X ) VCA;+l+2 X 

> 

(iii) for a > 1, 

ah'Kz) - fk{az) 



^ (q - l)c^+iafc+i/gfc 1 ^ p f ak+i^k^k , 1. ^ R /^ Ofc+iARfc , It 
[Ck+i + az){ack+i + az) Uk \ack+i + az ax J \Ck+i+az x 

(a - l)ck+iak+il3k 1^ „ / afc+i/SfcRfe 1 ^ A V"fc+iAR-fc It A 

H trKfc ^ \ In, ^ h -1„, 

ax Uk \ack+i + az ax J \Ck+i + az x / 

^ (« - l)Cfc+i«i+i/^fc (a - l)cfc+iafe+i^fc 



> 

where = R^ (^Pk-i)- properties are independent of R^ and therefore hold for n — )• oo. 

Assume now k = 1. For any sequence of bounded non-negative real numbers Zn, we have by ( |44l ) and 
the continuous mapping theorem IIBTTI . 

/l(^n)-/l(^n) ^^0. (56) 
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Thus, properties (i) — {in) of fi{z) also hold for fi{z). By Definition [2] and Theoremjs] these properties 
imply the uniqueness of positive solutions to the fixed point equations z = fi{z) and y = /i(y), and 



hence the uniqueness of solutions to (52 1 for A; = 1. Moreover, note that 



\fk{a) - fk{b)\ < ^I±l^||Rfc f\a - b\. (57) 

Hence, 

|ei (x,3i) - ei = |/i (ei - /i (ei {x,^i))\ 

< |/i(ei {x,p,))-fi{ei (x,^i))| + |/i(ei (x, ^i)) - /i (ei {x,p,))\ 
<en + ^^^\\Ki{Po)f\ei{xJ^)-ei{x,P^)\ (58) 
for some sequence of real numbers e„, satisfying e„ > 0. By Lemma 6 llRi (/3o)|| < M, almost 

n— 5>oo |_[ ^ ' 

surely, for some M > 0. Thus, for x < ^2^\^!^2 ^ri^ some 5 > 0, we have 

|ei(x,^i) -ei(x,^i)| < ^ ^^0. (59) 

Since ei {x,l3i) and ei (x,/3]^) are (almost surely) bounded on any closed subset of IR+ \ {0} and have 
analytic continuations for x G C\]R^, we have by Vitali's convergence theorem [I49II that the convergence 
holds for any X G ]R+ \ {0}. 

The last convergence implies by the continuous mapping theorem that, 

mi(x, (3i) — nil (x, /3i 0. (60) 

We now assume k = 2. The last convergence implies f2{z) — )■ f2{z), almost surely. The same steps can 
therefore be applied to show that m2{x,Pi) — m2 (x,/^^) > 0. This terminates the proof as this 
process can be iterated k times. 



n—^oo 



Appendix D 

Proof of Theorem |4) Fundamental equations 

The proof follows essentially the same steps as the proof of Theorem 2 in [7 1 and will not be given 
in full detail here. In order to prove the uniqueness of solutions {gk, 9k, ^k), it is sufficient to show by 
Theorem [S] that the A'-variate function h : (xi, . . . , xk) ^ {hi, . . . , hx) as defined below, is a standard 
interference function (see Definition |2]). For k = I, . . . , K, we define 

1 Sk Sk 

/ifc(xi,...,xx) ^ — V— 4^^^-^ (61) 
nk 9kSk,j0k 
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where 

5A: = ;^trT|QfcT|(xfcT|QfcT|+I„,) ' (62) 
and 6k, k = 1, . . . , K, form the unique jointly positive solution to the following fixed-point equations 

4 = ^trRjf;^^R, + ll^') \ (63) 

Nk Ok p J 

The only difference to ||7j| is the definition of gk- In our case, gk is directly defined as a function of 
Xk, whereas bk in fj] (using their notations) is given as the solution of another fixed point equation. 
However, the behavior of bk and gk as seen as functions of Xk is identical. In particular, let Xk > x'^ > 
and denote by gk and the corresponding values of ( |62] ), respectively. One can easily verify that the 
following conditions hold:(i) gk < g'k and {ii) Xkfjk > ^'kd'k- "^^'^ remaining steps are identical to [7| 
and will not be repeated here. By showing h(xi, . . . , xk) to be a standard interference function, we have 
proved by Theorem |8] that the following fixed-point algorithm, which iteratively computes 

4+' = hk{xf, . . . , 4^ k = l,...,K (64) 

for t > and some set of initial values xf\ . . . ,x^l^\ converges as t — ;> oo to the unique fixed point 
(51, • • ■,9k)- 

Appendix E 
Proof of Theorem [5j Mutual information 

The key idea is that the double-scattering model can be considered as the Kronecker channel model as 
considered in Theorem |9] with random correlation matrices. Assume now a Kronecker model, for which 
the matrices H;^. are given as 

Hfc = ^Z,,W2,fcT| (65) 



where e C^^^*" is a deterministic matrix and V^2,k and are defined as in (17). Further assume 
that limsup7v||Zfc|| < 00 for all k. Thus, we can apply Theorem [9] to obtain the following deterministic 
equivalent Vn{p) of In{p)' 

Vn{p) = - logdet + p^ekZkZ^ j^^N i^""' ^ efcT|QfcT|) - ^ J] ""^ekek 

k=l / k=l k=l 



(66) 
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where ek,ek,k = 1, . . . , K, are given as the unique solutions to the following equations 

efc = ^trT|Q,T|(efcT|Q,T|+I„,) ' 

efc = — tr Z,,ZH e^ZiZ^ + (67) 

such that efc, Cfc > for all k. Recall that the matrices are the covariance matrices of the channel inputs 
Xfc. Thus, the channel model is equivalent to a channel = -^^L^y^ 2,k^k' where = TsQ^, 
with uncorrelated channel inputs ~ CA/'(0, 1^^)- 

For the double-scattering channel, the matrices Z^ are random and defined as 

1 



^ VA^;^^^'''^^' ^^^^ 
Let (ri, T ^ P) be the probability space generating the random sequences of matrices (Wi ^(a;))^^^. There 

exists A C 0. with P{A) = 1, such that for each a; G A, we have limsup^||Zfc(a;)Zfc(cL!)'^|| < oo ( [1461 ). 
Thus, for each of these u, the matrices Z^Z^ satisfy the criteria of the correlation matrices of Theorem |9] 
Let {Q',T',P') the probability space generating the matrices 'W2^k- Thus, for every lo £ A, there exist a 
A'{u}) C Q' with P'{A') = 1, such that for all u' G A'{uj), Vn{p) is a deterministic equivalent of In{p)- 
Denote by (17 x x F',Q) the product space generating the matrices Wi fc(a;) and Wi fc(a;') and 

denote hy B <Z Q x ^' the space of all tuples {ui,ui'), such that u) ^ A and G A'{ijj). By the Fubini 
theorem, we have Q{B) = 1, which proves that Vn{p) — In{p) — ^ 0, almost surely. However, Vn{p) 
is a random quantity, which depends on the matrices Z^. Therefore, we will need to obtain an iterative 
deterministic equivalent In{p) of Vn{p)- 

The first step is to replace the fixed-point equations ([67]) that depend on Z^ by deterministic ones. Let 
us define the quantities efc,jj, efc,ij, for i G {1, . . . , K}, j G {1, . . . , N^}, which are given as the unique 
solutions to the following set of fixed-point equations: 

efc,ij = — trZfc^jjZ^jj ^ e^^j.jZ^^jjZyj J + -Iat (69) 

where 

'■^k,i,j 

• • • Zfcj_iZfcj_|_i . . . Zfc^TVfc] ; i = k 



Obviously, ek,i,j and ek,i,j are independent of the vector Zj j. In addition, we define 



Z = maxlimsup||ZfcZ^ ||, T = maxlimsup||Tfc| 



N 

n = mmn/j, c = — 
k n 



and 
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Oij = max \ek,i,j - ek\ , aij = max \ek,i,j - ek\ 



Thus, we have for N large. 



(70) 



Since the right-hand side (RHS) of the last inequality is independent of k, we have 



On the other hand, for i k, we have for N sufficiently large, 

I ^k,i,j 6fc I 



(71) 



— trZfcZ^ I ^ e^^jjZ^^jjZ^jj + -1^ 1 I ^(e^ — ee^ij)ZiZ^ + ei^ijZijz\^j | I ^e^Z^Z^ + -1^ 



n 



-1 



-1 



Zj^j I ^ e^Z^Zy + -Iat j ZfcZ^ I ^ e^ijZ^^jjZ^jj + -Iat j Zjj 



^Z^T 



n 



(72) 



where the last inequality is due to ek^ij < T and 
Similarly, one can show that 



< llzjjlpll A|| < ^||A||, for any matrix A. 



2 2- P^^^T pZ 

\ek,k,j - efcl = cAp z ttkj H ^ • 



n n 



It follows from (72 1, (73 1 and (71 1, that 



9r.9 P^Z'^T pZ 9 ^9^9 p^Z'^T pZ 

o-i j < cKp^Z^ai i + + — < cKp^Z^T^ai j + + — . 

n n ' n n 



(73) 



(74) 



Now, for any p < a/ cKZ^t^ ^^'^ e > 0, we have 



ai,j < — (1 + pZT) , ajj < 
en 



pZT^ 



en 



{1 + pZT) 
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Let /i = max{2^ (1 + pZT) , (1 + pZT)}, then we finally have 



(75) 

n n 

The last result establishes that for sufficiently small p, the differences between the solutions [e^^ij, ek,i,j) 



to ([69]) and the solutions (cfc, e^) to ( [67] ) are uniformly bounded by p and vanish as n — cxd. Moreover, 
efc and Ck^ij have an analytic continuation on z = ^ C \ and are uniformly bounded on all 
closed subsets of z G C \ Thus, in particular for p G IR+ \ {0} and all k, we have by the Vitali 
convergence theorem ll49l 







and hence 



Af-s>oo 



0. 



(76) 



As a consequence of (76 1, we can now write 



(a) J_ 
^ "-fc 

(c) 1 



'^k,j 



Zfcj 



E 



ftrRfc(E£ieu,fcZiZH + il 



H , li 



Af 



(d) 1 



E 



trRfc(E£ie.Z.ZH + il 



H , li 



N 



-1 



il + ^trR, Er=ie-.Z,ZH + ll 



AT 



(77) 



where (a) follows from (|76]l since 

-1 



\i=i P ) \%=\ ^ 



-1 

N I Zfcj 



< 



n 



Af-s>oo 



-^0, (78) 



(6) is due to Lemma|3j (c) is a consequence of Lemmas [2] and |4] and (d) is obtained by applying (|76]l a sec- 
ond time. Next, we would like to find deterministic equivalents of the terms ;^tr {Yld=i CjZjZ,^ + p^-Af 
We cannot directly apply Theorem |9] at this point since the are defined as functions of Z,fc. However, 
based on the relations ( [75] ) and ( [78] ), Theorem [9] (see llT4l Theorem 1]) can be shown to hold also for 
the matrix model under study. Thus, 



-^tr Rfc f V e,Z,zH + x A 



(79) 
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where fk tor k £ {1, . . . ,K} are defined as the unique solution to the following fixed-point equations 

JVfc ^ 1 + ekSkjJk 

fk = T^tr Rfc f V fiKi + -I^) (81) 
such that fk>0 for all k. Replacing ( [SO) ) in ( [8T] ) leads to 

A = ^tr R, f E TT^^ + -1^1 • (82) 



Thus, 



1 Skjfk 

Cfc = — > + efc 

nfc 1 + SkjekJk 



(83) 



where ek is a sequence of random variables, satisfying ek — — -)• 0. Consider now the following system 
of equations 

Sk = — trTfc fcfcTfc + In,^ 
Uk V 



Cfc = — > — ^-r + ek 

nk^l + SkjekJk 

/. = -^tr R, f R. X; + -In 



and its deterministic counterpart 

Qk = — trTfc (gk'ik + Infc^ 
rik 

SkjOk 

9k- ^ 



-t- 

rik ^ 1 



Define the quantities: 



1 , K _ ^ \ ^ / K _ 



-1 



71 = max|efc - , 72 = max |efc - fiffc| , 73 = max - 4| , e = max|efc|. 

K k k k 

Straight-forward calculations lead to the following bounds: 

71 < c573 + cS'^R^p^-12 + e, 72 < 71 ^^ 73 < KSR^p^^2 + KS^R^T^p'^ji. 
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Combining these results and using the fact that e — '—^ yields for p sufficiently small, 

N^oo 



71,72,73 



-> 0. 



(84) 



Since ek,gk,ek,gk, fk^^k are all (almost surely) bounded for p in any closed subset of \ {0} and 
have analytic continuations for p G C\]R^, we have by Vitali's convergence theorem P9]| that ( |84| ) holds 

for any p € M+ \ {0}. 



Coming now back to Vn{p) as given in (66 1, we have from the continuous mapping theorem |31] that 

K 



k=l 



logdet I„, + efcT/ 



rikekek 



Moreover, since |Efc=i (efc - Qk) ZfcZ^ 



log det [In^ + gkTk ) - UkQkQk 
0, we have 



7V-s>oo 



0. 



(85) 



^logdet ^Ijv + /o J^efcZfcZ^^ - ^logdet ^Ijv + p 5fcZfcZ^^ 
Applying Theorem [9] to the last term yields 
^ log det (^n + pY,^ 5fcZfcZ^^ 



-> 0. 



(86) 



K 



log det f Iat + p ^ ^ ^^Rfe j - V7 XI log det (Iat, + gkSkSk) + UkQkQk > 0. (87) 



fc=i 



Nk h 



N 



k=l 



Combining ( 85 1 and ([87]) finally leads to 



Vn{p)-In{p)-^0 

N^oo 



(88) 



where 



K 



In{p) = ^logdet ( In + p^^^^Y^Kk 1 + ^ ^ log det (Iat, + gkSkSk) - UkQkQk 



K 



N 



k=l 



Nk Sk 



N 



k=l 



1 ^ 

k=l 



nk9k9k- 



(89) 



This concludes the proof of part {i). 

In order to show the convergence in the mean (part (ii)), we will pursue the same approach as in IH, 
Q. Define the following functions: 



K 



-1 



K 



-1 



m;v(^) = ^tr J]H,hH-zI 

\fc=i 

One can show that 

E[lNip)]-lNip)-- 



N 



- ^ \ . ( "fc 9k9k r, T 



N 



\k=l 



Nk Sk 



N 



E[mjv(-u;)] 



mN{-oj) 



doj. 



(90) 
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Since both m]^{uj) and miv(w) are uniformly bounded by ^, it follows from dominated convergence 
arguments, Theorem |9] and ([84]) that, for all a; > 0, 



E[miv(-6^)] 



UJ 



mN[-uj) 



(91) 



Moreover, 
1 



E[mN{-uj)] 



UJ 



mN[-uj) 



< 



1 

UJ 

1 / 1 



E[mN{-uj)] 

K 



+ 



UJ 



mN[-uj) 



.k=l 



K 



nk gk9k 



< 



2KRST 



(92) 



where = max^ limsup||Rfc||, S = max^ limsup||Sfc||, T = max^ limsup||TfcQfc||. Since 
finite and integrable over oo), it follows from the dominated convergence theorem that 

E[In{p)]-In{p)- ^0. 

N^oo 

Appendix F 

Proof of Theorem[6J Optimal power allocation 

The proof follows closely those of |[25l Proposition 5] and llT4l Proposition 3]. 
We first recall the following property of concave functions (see e.g. lISOll ): 



2KRST 



IS 



(93) 



Property 1: A function / (Qi, . . . , Qi^) is strictly concave in the Hermitian nonnegative matrices 
Qi, . . . ,Qk, if and only if, for any couples (Qi^, QiJ , . . . , (Qa'„, QirJ of Hermitian nonnegative 
matrices, the function 

0(A) = / (AQi„ + (1 - A)Qib, . . . , AQk„ + (1 - X)QKb) , A G [0, 1] 

is strictly concave. 

Consider now In{p) seen as a function of A for Qfc = XQ,k^ ~ (1 ~ X)Q,k,,, where Qfc^jQfc^ are 
Hermitian nonnegative definite matrices, for A; = 1, . . . , i^. Thus, by the chain rule of differentiation. 



dlN{p) _ dlwip) ^ dlNjp) dgk ^ dlNjp) dgk ^ dlNjp) 05 
^ dgk d\ dgk dX d5k dX' 



dX 



dX 



k=l 



One can verify that the partial derivatives of {p) with respect to gk,gk,5k, respectively, satisfy 

OIn ip) diN (p) OIn {p) 



dgk 



dgk 



d5k 



0, 



(94) 



(95) 
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due to the defining relation {TA). Thus, 

K 



k=l 



-1 1 



— = ^ = E (I.. + ,.T| (AQ.. + (1 - A)QO T|) ^ T| (Q.^ - Q,, 



(96) 



The second derivative therefore reads 



dA2 



K 



k=l 



In, + 5fcT| (AQfc„ + (1 - A)QfcJ (Qfc„ - Q 



(97) 



where are Hermitian nonnegative definite and are Hermitian. Let = Uj^DfeU^ be the eigen- 
value decomposition of A^, where Ufc G ^i^kxnk unitary matrices and = diag . . . , d^^nk) ^ 
0. Moreover, denote = B^Ufc = [z^ i . . . Then, 

^tr [AfcBfe]2 = ^trDfcU^^BfeAfcBfeU = ^trDfcZ^AfcZfc = ^ ^ 4,iZfcj AfcZ^^ > 0. (98) 



N 



N 



N 



If Tfc ^ 0, or equivalently, if is invertible, we have A^ >- and ( [98] ) holds with strict inequality for 
Qfc„ 7^ Qfcf Thus, if any of the matrices is invertible and Q^^ ^ Q^^,, we have 



^ 1 

j;-tr [A,Bfc]^>0 



(99) 



fe=i 



and hence ^ ^"i^-* < 0. Thus, In{p) is strictly concave in the matrices Q^. Due to (95 1 it is then sufficient 



to maximize 



(100) 



logdet +fffcT^QfcT = 

with respect to and with the constraint — trQ^ < P^. The solution to this problem is well-known 
Bn and given by the water-filling solution stated in the theorem. However, since depends on Q/^., 
this solution must be computed iteratively by Algorithm 1. 



Appendix G 

Proof of Theorem[7J SINR of the MMSE detector 

Similar to the proof of Theorem [5] let us consider the matrix model 

1 



where 



Z, 



(101) 



(102) 
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Then, 



N 

7k,j = Pk 



H Hi ^ H H i 

,jtk,j ^2k]^k / HjHj — PkjtkJ ZfcW2,fcjW2,fejZ;j + -ItV ZfcW2,fco-. (103) 

nk nk P J 

It was shown in the proof of Theorem [s] that, almost surely, limsup;Y||ZfcZ^ || < 0. From the Fubini 
theorem. Lemma [2] and Lemma |4) we therefore have 



N 

lk,j - Pk 



jtkj — trZfeZ^ ^HjH,^ +/'Infe 



0. 



Applying Theorem |9] (i) leads to 



N 

lk,j - Pk,j 



tfcj^trZfcZ^ (^^eiZiZ^ + pln}j 



-1 



Af-5-oo 



-> 



where Cj are given as the unique solutions to ( [67] ). Notice now from ([67]) that 

-1 



and that max^ | — i^a: I 



N^oo 



ek = — trZfcZ^ fVeiZiZ,^ +/9l„,^ 
> by (l84]). This finally implies 



M . a.s „ 

7fe,j - Pk,jtk,j9k > 0. 



(104) 



(105) 



(106) 



(107) 



Appendix H 

Proof of Corollary[2J Sum-rate with MMSE decoding 
Part (i) is a simple consequence of Theorem [t] and the continuous mapping theorem. 
For Part {ii), first notice that Rn{p) < In{p) and Rn{p) < In{p)- Thus, 



K n 



^n{p) - j^Yl^ log (1 + Pkjtkjgk] 



k=i j=i 



< I Nip) + Iwip) = (pN^p)- 



Since E [4>n{p)] < oo by Theorem [S] (ii), it follows from dominated convergence arguments that 

E[RNip)]-RNip) ^0. 

Appendix I 

Proof of Corollary O Rayleigh product channel 

Under the assumptions of the corollary, the fundamental equations in Theorem |4] reduce to 

1 

9 

9 
5 



(108) 



(109) 



1+5 
S 6 
Nl + g6 
1 

^ s ~^ N p 



(110) 
(111) 
(112) 
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From (fTTOll, we have 



Solving for 6 and replacing g by ( |113| ) yields 



1) 

Solving ( |112| ) for (5 and replacing 5 by ( |113| ) leads to 

l-i^(l-g) 
= • 

N p 



Equating ( 114i and (115 1 and rearranging the terms as a polynomial in g finally yields 



9 -9 



S_ 
N 



1 

K 



+ 5 1 



S_ 
N 



1 

K 



S 



1 + 



1 



S 1 



NK p 



(113) 



(114) 



(115) 



(116) 



By Theorem [4| only one of the roots of this polynomial satisfies g,g,6 > 0. Now, ( |1 13| ) implies g < 1, 
( |114| ) implies 5 > 1 - ^ ( |115| ) implies ^ > 1 - ^. Hence 5 G (l - min [;^, ;f ] , l). 
Similarly, /jv {p) reduces under the assumptions of the corollary to 



In (p) = log 1 + yo 



NKgg\ , KS 



+ — log (1 +g6) + K log (1 + <?) - 2Kgg. 



(117) 



Replacing | by (7 + — 1 in the first term, 5 by ( |1 14[ ) in the second term, g by ( |113| ) in the third term 
and 59 by (1 — ^) in the last term leads to the desired result. 



The simplification of Theorem |7| is immediate since 7^ = Phji^jg^ 



by ( fTT3| ). 



Appendix J 

Matlab code related to the multihop AF MIMO relay channel 



A. Code for Corollary^ corollaryl 



m 



function I_k = corollaryl (k, K, a, rho, c) 

% Compute deterministic equivalent I_k of the asymptotic mutual information 
% Input parameters : 

denotes which function I_k to compute 
total number number of hops 
vector containing the path loss factors 
vector containing the power budgets rho_k 
vector containing the matrix dimension ratios 



% k 
% K 

% a = [alpha_l , . . . , alpha_{ K} ] 
% rho = [rho_0, . . . , rho_{K-l } ] 
% c = [c_l, . . . , c_K] 
% Calculate asymptotic power normalization factors (Lemma 1) 
b = zeros ( 1 , k) ; 
b (1) = rho (1) ; 
for 1=2 : k 

b(i) = rho(i) / ( 1+a ( i-1 ) *rho (i-1 ) ) ; 

end 

% Calculate capacity 

I_k = 1 /K* (theorem2 (k,l,a,b,c) - theorem2 (k, l,a, [0,b(2:end) ] ,c) ) ; 
end 



B. Code for Theorem^ theorem2 . 



m 



denotes which function J_k to compute 
argument of J_k 

vector containing the path loss factors 
power normalization factors 

vector containing the matrix dimension ratios 



function J = theorem2 (k, x, a, b, c) 

% Recursively computes deterministic equivalent of J_k 
% Input parameters : 
% k 
% X 

% a = [ alpha_l , . . . , alpha_{ k } ] 
% b = [beta_0, . . . ,beta_{k-l} ] 
% c = [ c_l , . . . , c_K] 
if (k==l) % J_l is given in closed form 
["■,e] = theorems (k-1, X, a, b, c) ; 

J = c(k)*log(l + x*a (k) *b (k) / (c (k) +e) ) + log (1+e/c (k) ) - e/(c(k)+e); 

else % J_k, k>l, must be computed recursively 

["■,e] = theorems (k-1, X, a, b, c) ; 

J = c (k) *theorem2 (k-1, x*a (k) *b (k) / (c (k) +x*a (k) *b (k) + e),a,b,c)... 

+ c(k)*log(l + x*a (k) *b (k) / (c (k) +e) ) + log(l + e/c (k) ) - e/(c(k)+e) 

end 
end 
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C. Code for Theorem^ theorems. m 

function [m,e] = theorems (k, x, a, b, c) 

% Recursively computes the quantities e_k 
% Input parameters : 

% k : denotes which function e_k to compute 

% X : argument of e_k 

% a = [alpha_l, . . .,alpha_{k}] : vector containing the path loss factors 
% b = [beta_0 , . . . , beta_{ k } ] : power normalization factors 

% c = [c_l, . . . , c_K] : vector containing the matrix dimension ratios 

if (k==0) % For k=0, e_o and m_0 are given in closed form 

e = - (x*a (1) *b (1) * (1-c (1) ) +c (1) ) /2 + sqrt ( (x*a ( 1 ) *b ( 1 ) * ( 1-c ( 1 ) ) +c ( 1 ) ) " 2 . . . 
+ 4*x*a (1) *b (1) *c (1) "2) /2; 

m=c(l) / (a (1) *b (1) / (c (1) +e) +l/x) + (l-c(l))*x; 
else % For k>0, e_k is given as the solution to a fixed point equation 

e = 0; 

eold = 1; 

while abs (e-eold) >le-6 % the error tolerance le-6 can be changed 
eold = e; 

e = c (k+1) * (c (k+1) +eold) - c (k + 1 ) * ( c (k+1 ) +eold) ~ 2 / (x*a (k+1 ) *b (k+1 ) ) . . . 
* theorems (k-1 , x*a (k+1 ) *b (k+1 ) / (c (k+1 ) +x*a (k+1 ) *b (k+1 ) +eold) , a, b, c) ; 

end 

m = x*c (k+1) / (c (k+1) +e) ; 

end 
end 
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